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ABSTRACT 

Kashlinsky et al. (2005) find a significant cosmic infrared background fluctuation ex¬ 
cess on angular scales > 50 arcsec that cannot be explained by instrumental noise or 
local foregrounds. The excess has been tentatively attributed to emission from pri¬ 
mordial very massive (PopIII) stars formed < 200 Myr after the Big Bang. Using 
an evolutionary model motivated by independent observations and including various 
feedback processes, we find that PopIII stars can contribute < 40% of the total back¬ 
ground intensity ~ I — 2 nW m~^ sr“^ in the 0.8-8 p,m range) produced by 
all galaxies (hosting both PopIII and PopII stars) at 2 ; > 5. The infrared fluctuation 
excess is instead very precisely accounted by the clustering signal of galaxies at z > 5, 
predominantly hosting PopII stars with masses and properties similar to the present 
ones. 
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1 INTRODUCTION 

Observations of the infrared background provide important 
information on the emission of cosmic Inminous sources 
throughout the history of the Universe. It has been sug¬ 
gested (Santos, Bromm & Kamionkowski 2003; Salvaterra 
& Ferrara 2003) that a large fraction of the measured Near- 
InfraRed (1-10 /rm) cosmic Background (NIRB) arises from 
redshifted Lya line photons and nebular emission produced 
by the first very massive metal-free stars. This hypothesis, 
however, is very demanding in terms of the required conver¬ 
sion efficiency of baryons into stars (Madau & Silk 2005). A 
large NIRB contribution from such stars has more recently 
been rejected by the paucity (< 3) of 2 ~ 10 candidate 
sources in Hubble Space Telescope ultra-deep observations 
(Salvaterra & Ferrara 2005). Nevertheless, a more modest 
contribution from very high redshift galaxies, whose clus¬ 
tering should leave a distinct signature on small-scale an¬ 
gular fluctuations of the background light (Magliocchetti, 
Salvaterra & Ferrara 2003; Kashlinsky et al. 2004; Cooray 
et al. 2004), is still possible. 

Kashlinsky et al. (2005) have recently found signifi¬ 
cant NIRB fluctuations in deep exposure data obtained with 
Spitzer/IRAC (Fazio et al. 2004a, 2004b) in four channels 
(3.6, 4.5, 5.8, and 8 /rm), after Galactic stars and galaxies 
bright enough to be individually resolved by the instrument 
have been carefully subtracted. With the only exception of 
the 8 ^m channel, the shape and amplitude of the power 


spectrum cannot be reproduced by either contributions from 
intervening dusty. Galactic neutral hydrogen gas (cirrus) 
or from local interplanetary dust (zodiacal light). Ordinary 
galaxies (z < 5) produce fluctuations due to their clustering 
and shot-noise. The faint flux limits (> 0.3 /rJy) of Spitzer 
data allow to push their residual clustering contribution be¬ 
low the level of the excess signal at relatively large (^50 arc- 
sec) angular scales (Kashlinsky et al. 2005). The shot noise 
component, estimated directly from galaxy counts, fits the 
observed fluctuations at smaller angular scales, and rapidly 
fades away at larger angles. The residual large scale signal 
has been ascribed by Kashlinsky et al. (2005) as coming 
from very distant (2 > 5) sources provided their total flux 
contribution is > 1 nW sr“^. The aim of this Letter is 
to show that this is indeed the case. 

The layout of the paper is as follows: in Section 2 we 
will briefly describe the adopted model, while in Section 3 we 
provide predictions for the NIRB intensity and fluctuations 
and compare the latter ones with the results of Kashlinsky 
et al. (2005). Section 4 summarizes our conclusions. 


2 THE MODEL 

Schneider et al. (2005) have presented a scenario for the 
formation of galaxies in a concordance AGDM cosmolog- 
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ical model* which includes a self-consistent treatment of 
two key feedback processes: (i) radiative feedback, suppress¬ 
ing star formation in H 2 -cooling halos and the formation of 
low-mass galaxies due to the effects of UV background radi¬ 
ation approaching the reionization epoch, and (ii) chemical 
feedback, which controls the transition from metal-free stars 
(PopIII) to ordinary stars (PopII) through the progressive 
enrichment of star forming gas with heavy elements released 
by supernova explosions (Schneider et al. 2002, 2004, 2005; 
Bromm et al. 2001). Chemical feedback propagates through 
the hierarchy of galaxy mergers from progenitors to their 
descendants so that, at each redshift, existing halos which 
are allowed to form stars are classified as PopII (PopIII) 
galaxies depending on whether the halo itself or any of its 
progenitors have (have not) already experienced an episode 
of star formation. 

Within this model we can compute the comoving spe¬ 
cific emissivity, e,^, which is given by 

r°° l-MmaxW) .2 

e^{z) = / dzl,,{txy) / M„——-—{Mh,z')dMh,{l) 

Jz J aMhaz 

where d^n/dMudz is the formation rate of halos of total 
(dark+baryonic) mass Mh with corresponding stellar mass 
M*, and Mmin{z) is the minimum mass of halos (corre¬ 
sponding to a virial temperature of lO'^ K, i.e. Mmin{z) ~ 

10® Mq (1 -I- z/10)“®^^) allowed to form stars at redshift 
Z', Mmax{z) is the maximum mass of halos which depends 
on the details of their merging history (see Schneider et al. 
2005); lv{tx,z') is the template specific luminosity for a stel¬ 
lar population of age tz,z' (time elapsed between redshift z' 
and z). Following the results of Schneider et al. 2005, both 
PopII and PopIII stars are assumed to form according to 
a Salpeter Initial Mass Function (IMF) with an efficiency 
/* = 0.1. The emission properties of PopII stars are taken 
from the GALAXEV library (masses in the range 0.1-100 
Mq, metallicity Z = 10“^^©; Bruzual & Chariot 2003) and 
those of PopIII stars are based on the synthesis model for 
metal-free stars of Schaerer (1-100 Mq, Z = Q-, Schaerer 
2002 ). 

A set of independent observational constraints can be 
accommodated within this model (Schneider et al. 2005): 
the number counts of dropout galaxies at a ~ 6 (i-dropouts; 
Bouwens et al. 2005b) and at z ~ 10 (J-dropouts; Bouwens 
et al. 2005a), and the value of the optical depth for Thomson 
scattering = 0.16±0.04 measured by the WMAP satellite 
(Kogut et al. 2003). Indeed, the predicted surface density of 
i-dropouts at a limiting magnitude of i = 28 is ~ 5 arcmin”^ 
(observed value ~ 4.7 arcmin“^; Bouwens et al. 2005b). The 
model is also consistent with the three candidate dropouts 
at 2 ~ 10 reported by Bouwens et al. (2005a). Moreover, the 
large value of Te can be matched with standard values of the 
ionizing photon escape fraction from galaxies (/esc < 0.2) 
and gas clumping factor C ~ 10, resulting in a reionization 
redshift of 13.2 (Schneider et al. 2005). 


* We adopt a ACDM cosmological model with parameters = 
0.3, Oa = 0-7, h = 0.7, fig = 0.04, n = 1, and erg = 0.9, which lies 
within the experimental errorbars of WMAP experiment (Spergel 
et al. 2003) 



Figure 1. IR background from sources at 2 > 5 in the wavelength 
range 0.3-8 pm. The dotted (dashed) line shows the contribution 
of PopII (PopIII) galaxies. The sum of the two components is 
denoted by the solid line. The sharp drop below 0.8 pm is due to 
the absorption of the intergalactic medium for A < Anya in the 
emitter rest frame. 


3 NIRB INTENSITY AND FLUCTUATIONS 

The model presented in Section 2 (referred to as model C in 
Schneider et al. 2005) allows to compute the NIRB intensity 
Ji/Q from 2 > 5 galaxies seen at frequency uo by an observer 
at redshift 2 o = 0. This can be done in a fashion similar to 
Salvaterra & Ferrara (2003) as: 

■^- 0 (^ 0 ) = / eAz)e-^^ff^''°-^°’^^^dz, ( 2 ) 

J ZQ 

where u = vo{l ± 2 )/(l ± 20 ), dl/dz is the proper line el¬ 
ement, re//(j^o, 2 o, 2 ) is the effective optical depth at vq of 
the intergalactic medium between redshift 20 and 2 (see Sec¬ 
tion 2.2 of Salvaterra & Ferrara 2003 for a full description of 
the IGM modelling), and £ 1 ^( 2 ) is provided by equation (1). 
Our results are shown in Figure 1. The background intensity, 
vJv ~ 1 — 2 nW m~^ sr“^, is almost constant in the 0.8-8 pm 
range. PopII galaxies dominate the NIRB in the entire wave¬ 
length range, while PopIII galaxies contribute at most 40% 
of the total intensity (at A 1.5 pm), via their strong Lya 
line emission. We have checked that the sources responsible 
for such emission are too faint to be resolved and identified 
at the Spitzer magnitude limit in all four channels. The total 
contribution of all galaxies at 2 > 5 is 1.55, 1.45, 1.07 and 
0.74 nW sr“^ in the 3.6, 4.5, 5.8 and 8.0 IRAG bands, 
respectively, representing a substantial fraction (~ 10-20%) 
of the total NIRB intensity estimated from integration of 
Spitzer galaxy number counts (Fazio et al. 2004b). 

The angular correlation function of intensity fluctua¬ 
tions 5J due to inhomogeneities in the space distribution of 
unresolved sources (i.e. with fluxes fainter than some thresh- 
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old Sd) is defined as: 

C{e) = {5j{e',cj>')5j{e'',cj>'')), ( 3 ) 

where {O', 4>') and {9", (j>") identify two positions on the sky 
separated by an angle 9. The above expression can be writ¬ 
ten as the sum of two terms, Cp and Cc, the first one 
due to Poisson noise (i.e. fluctuations given by randomly 
distributed objects), and the second one owing to source 
clustering. It can be shown that the shot noise contribution 
originating from z > 5 galaxies is negligible, so we only con¬ 
centrate on fluctuations which stem from the clustering of 
these sources, i.e. we assume C{9) = Cc- 

The method adopted here is similar to that presented 
in Magliocchetti et al. (2003), whereby angular fluctuations 
are obtained by means of the expression: 



{i + zF 





X 


du ^g{r,z) 


( 4 ) 



10 100 10 100 
27 T/q [arcsec] 


where x = 1{1 + z) is the comoving coordinate, r = -|- 

^ 292 -^ 1/2 universe and in the small angle approxi¬ 

mation), and tv is defined as in Section 2. 

The spatial two-point correlation function of a class of 
galaxies ’g\ ^g{r,z), in general results from a complicated 
interplay between the clustering properties of the underly¬ 
ing dark matter and physical processes associated to the 
formation of such galaxies (see e.g. Magliocchetti & Por- 
ciani 2003). However, the sources we are considering in this 
work are small enough (typical dark matter masses around 
10® — 10 ^®Mq for z ~ 5, values which rapidly decrease when 
moving to higher redshifts, see Schneider et al. 2005 for 
further details) and set at high enough redshifts to ensure 
that - if present - their sub-halo behavior falls into angular 
scales which are too small to be detected. We can there¬ 
fore safely assume a one-to-one correspondence between ha¬ 
los and galaxies and write: ^g{r,z) = ^(r, z)bgg(z), where 
5(r, 2 ) is the mass-mass correlation function (evaluated by 
following the Peackock & Dodds 1996 approach up to 100 
Mpc comoving) and beg(z) the bias associated to all dark 
matter halos massive enough to host a galaxy at redshift 2 . 
This latter quantity was then obtained respectively for PopII 
and PopIII galaxies by integrating the function b{Mh,z) - 
representing the bias of individual haloes of mass Mh - op¬ 
portunely weighted by the number density of such sources 
dn/dMh{Mh,z) f{Mh,z) as: 


rj T\/f / ■> -r ’ 


where dn/dMh{Mh,z) is the Sheth & Tormen (1999) mass 
function of the dark matter halos, f{Mh,z) is the fraction 
of halos that at any mass and redshift host a PopII (PopIII) 
galaxy and Mmin{z), Mmax{z) are the limiting mass values 
for such halos to host either a PopII or a PopIII source (see 
Section 2). 

The above prescriptions then allow to estimate C{9), 
i.e. the contribution from the clustering of unresolved 2 > 5 
(PopII and PopIII) galaxies to the background fluctuations 
at different wavelengths. An important feature to mention 


Figure 2. Power spectrum of fluctuations at the IRAC bands. 
Different symbols show data of different observed fields (Kashlin- 
sky et al. 2005): filled circles corresponds to the QSO 1700 field, 
open circles to EGS. Data are averaged over wide bins to increase 

_/2 

the signal-to-noise ratio. The errors are ' , where Nq is the 
number of Fourier elements at the given g-bin. Dotted (dashed) 
lines show the contribution from PopII (PopIII) galaxies. Dot- 
dashed lines show the contribution of shot noise from remaining 
galaxies fainter than the limiting magnitude (Kashlinsky et al. 
2005). The solid line is the sum of the different components. 


about the resulting C{9) is that - independent of the con¬ 
sidered frequency and only due to the high redshift of the 
sources in exam - it features a sharp drop at around 300 
arcsecond beyond which it rapidly approaches zero. This im¬ 
plies that signals on scales larger than the above value such 
as those detected by some experiments (e.g. Hauser & Dwek 
2001; Matsumoto et al. 2005) have to rely on the presence 
of a more local population of unresolved sources. 

Finally, we can move to the calculation of the power 
spectrum P{q) originating from the above fluctuations, 

/»oo 

P{q) = 2-K C{9)Jo{9,q)9d9, (6) 

Jo 

where Jo is the zero-th order Bessel function. 

In order to confront our predictions with the Kashlinsky 
et al. (2005) results, P{q) has been evaluated in the 3.6, 4.5, 
5.8 and 8 fim IRAC bands. Our findings are summarized in 
Figure 2, where the dotted lines represent the fluctuations 
originating from PopII galaxies, the dashed lines those due 
to PopIII sources and the dot-dashed curves indicate the 
shot noise contribution from galaxies fainter than the lim¬ 
iting magnitude as estimated by Kashlinsky et al. (2005). 
It is clear from the plot that at the lowest frequencies the 
contribution of 2 > 5 galaxies can provide an excellent fit 
to the observed fluctuations; in particular it dominates over 
the shot-noise component at all angular scales greater than 
~ 50 arcsec. At 5.8 /rm, their importance is reduced but the 
presence of 2 > 5 sources is still required to properly account 
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for the amplitude of the signal at the largest probed scales. 
Finally, high redshift galaxies and shot-noise alone cannot 
reproduce the power spectrum at 8 /rm, where the cirrus 
probably becomes the dominant component (Kashlinsky et 
al. 2005). 

As a last point it is worth mentioning that, as expected 
within the present scenario, the contribution of PopIII galax¬ 
ies is negligible in all IRAC bands; the amplitude of their 
signal is in fact ^ 50 times lower than that of PopII galaxies 
in the most favorable case (3.6 pm), while the PopIII-to- 
PopII contribution ratio can go as low as 10“® at the highest 
frequency probed by IRAC. 


4 CONCLUSIONS 

Using a physically-motivated, observationally-tested model 
of the early Universe (Schneider et al. 2005), we compute 
the expected background radiation in the NIR by sources 
forming when the Universe was < 1 Gyr old. We find that 
the background intensity, ^1 — 2 nW m~^ sr“^, is 
almost constant in the 0.8-8 pm range. PopII galaxies domi¬ 
nate the NIRB in the entire wavelength range, while PopIII 
galaxies contribute at most 40% of the total intensity (at 
A ^ 1.5 pm), via their strong Lya line emission. Finally, we 
found that the infrared fluctuation excess on angular scales 
> 50 arcsec detected by Spitzer/IRAC (Kashlinsky et al. 
2005) is accounted very precisely by the clustering signal of 
galaxies at z > 5 predominantly hosting stars with masses 
and properties similar to the present ones. 

Two additional points are worth noticing: (i) very mas¬ 
sive stars (M > lOOM©) do not need to be invoked to ex¬ 
plain NIRB fluctuations and reionization history; (ii) be¬ 
cause of their small contribution {P{q) < 10“^® nW^ m“^ 
sr“^) to the observed power spectrum in all channels, ex¬ 
tracting the signal of the (very) hrst PopIII stars is ex¬ 
tremely challenging. Future instruments (as the James Webb 
Space Telescope) will be able to directly identify these 
sources up to « = 10 or above. Finally, the intensity of the 
NIRB provided by z > 5 galaxies falls short of accounting 
for the excess measured by IRTS (Matsumoto et al. 2005) 
and DIRBE (Hauser & Dwek 2001) experiments. The ori¬ 
gin of this component remains very puzzling (Salvaterra & 
Ferrara 2005) and might require either a revision of current 
model of zodiacal light subtraction or the existence of a large 
population of faint galaxies located at 2 ; ~ 2 — 3 (or both). 
Important insights on these issues are expected from the 
upcoming GIBER experiment (Bock et al. 2005), that will 
be able to simultaneously measure the total NIRB inten¬ 
sity and fluctuation power spectrum in the poorly known 
wavelength range 0.8-2 /rm. Such instrument, in addition, 
will allow a clear separation of the cosmological signal from 
local foregrounds (i.e. zodiacal light). 
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